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Abstract 

\ Similarity symmetries of the factorization chains for one-dimensional differential and 

Q> ■ finite-difference Schrodinger equations are discussed. Properties of the potentials defined 

by self-similar reductions of these chains are reviewed. In particular, their algebraic 
QhI structure, relations to (/-special functions, infinite soliton systems, supersymmetry, coher- 

' ent states, orthogonal polynomials, one-dimensional Ising chains and random matrices are 

(N 



o 



outlined. 
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1. Introduction. Spectral problems of various types play an important role in classical and 
psj ■ quantum mechanics and in the theory of integrable systems. Depending on the situation they 

O , may be based upon differential, finite-difference or integral equations, etc. Roughly speaking all 

^ ! such problems can be approximated by a search of eigenvalues of some matrices |l[]. This leads 

^\ \ to a number of common features in the methods of treatment of different spectral problems. 

We choose as a basic object of discussion the ordinary one- dimensional Schrodinger equation 



LiIj{x) = -i/j^^ix) -^u{x)iIj{x) = Xipix). (1) 



The eigenvalue problem for the operator L is completely defined by imposing some bound- 
^ ■ ary conditions upon iIj{x). However, it is convenient to consider (|I]) without such additional 

constraints and assume that x, A G C. When a non-formal meaning of the eigenvalues and 
^ . eigenfunctions will be necessary, we assume that ip{x) e L'^{R). 

There are several connections of (|l]) with finite-difference equations and discrete systems 
— the main subject of the present workshop. First, the scheme of treatment of (|l|) to be 
described below is universal, e.g., it can be easily reformulated for the second order difference 
equation. Second, finite-difference equations appear already in the consideration of (|l]) for 
different potentials. In the space of all equations of this form finite-differences enter through 
the symmetry transformations mapping Schrodinger equations onto each other: 

u{x) u{x), ip{x) — > ip{x), A ^ A. 

Third, in some cases the discrete spectrum eigenfunctions satisfy simple discrete equations. 
E.g., for u{x) oc one gets from the condition of square integrability of eigenfunctions A„ oc n, 
ipnix) OC Hn{x)e~^ where Hn{x) are Hermite polynomials satisfying three-term recurrence 
relation. Fourth, there is a direct relation of (|1]) to some lattice models of statistical mechanics. 

The paper is organized as follows. In the next several sections we review briefly some old 
results on self-similar potentials, their coherent states and the relation to g-special functions. 

^Invited talk by the second author at the Third International Workshop on Symmetries and Integrability of 
Difference Equations (16-22 May, 1998, Sabaudia, Italy). 
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Then we outline briefly a generalization of the construction to orthogonal polynomials. In the 
last two sections we present our recent results on the appearance of one-dimensional Ising chains 
and random matrices in the context of the Schrodinger equation (|I|). The list of references given 
in the end does not pretend to be complete, it contains mostly the papers which influenced 
somehow the work of authors. 

2. Self-similar spectra. Under the similarity transformations in (i- dimensional Euclidean 
space one counts usually the rotations, dilations and shifts of the coordinates, Xi QijXj + ai, 
i,j = l,...,d, i.e. the affince transformations. In the d = 1 case this is reduced to x —>■ qx + a, 
which will be used below under the assumption that the parameters q and a are fixed (for q 1 
one may set a = without loss of generality). 

Let us consider separately the effects of translations and scalings upon the Schrodinger 
operator L with potential u{x), eigenfunctions {ipj{x)} and the spectrum {Xj} (the use of the 
subscript j for notation of the spectrum is symbolic). After application of some symmetry 
transformation we get the operator L = —d'^/dx'^ + u{x) with the potential u{x), unnormalized 
eigenfunctions {ipjix)} and the eigenvalues {Xj}- So, the shift x ^ x + a brings in 

u{x) = u{x + a), {i)j{x)} = {ipj{x + a)}, {Xj} = {Xj}. 

The spectrum does not change under this transformation due to the absence of distinguished 
points in the boundary conditions. In the case of scaling transformations one replaces in (|l|) 
the variable x by qx, then multiplies the resulting equation by q^ in order to bring it to the 
canonical form and gets the operator L with 

u{x) = q^u{qx), {'tpj{x)} = {^{qx)}, {Xj} = {q^Xj}. 

Now one gets a scaled form of the spectrum A q^X. For |g| > 1 the spectral points are 
stretched with respect to the A = point, and, vice versa, for |g| < 1 we get a squeezed 
spectrum. The complete affine transformation of the spectral parameter requires joining of a 
shift X ^ X + h which is easily reached by the shift of the potential 

u{x)=u{x) + h, {'4'j{x)} = {^j{x)}, {Xj} = {Xj + h}. 

The combination of these three transformations gives 

u{x) = q'^u{qx + a) + h, {^jjj{x)} = {'ipj{qx + a)}, {Xj} = {q'^Xj + h}, 

i.e. 'trivial' similarity transformations of the spectral data. The key constructive idea, bringing 
in a non-trivial content into the above considerations, consists in a sticking to the systems with 
self-similar spectra for which {Xj} = {Xj}. Note that such a constraint is applicable to any 
spectral problem, not just to the Schrodinger equation. 

Analyze first the q = 1 case. The condition {Xj + h} = {Xj} shows that the discrete 
spectrum should appear in the form of a number of independent bilateral arithmetic progressions 
with the increment equal to h. In the simplest case the number of such progression is finite 
(= A^) and one can parametrize the eigenvalues as solutions of the equation Xj + h = A^+at or 
^pN+k = hp + Xk, k = 1, . . . N, p E Z. The continuous spectrum may appear in the form of 
arithmetic progressions of permitted bands. For instance, if one requires that the addition of h 
to potential is equivalent to an isospectral transformation x x + a, i.e. u{x) + h = u{x + a), 
then u{x) is a sum of a periodic potential Up{x + a) = Up{x) and the Airy part oc x known to 
lead to continuous spectrum filling the whole line, u{x) = Up{x) + hx/a. 
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Consider now the {q^Xi} = {Xi} case. Formally, the discrete spectrum appears in the form 
of bilateral geometric progressions accumulating near the A = point from below or above. 
Enumeration of pure point spectra in case of the finite number of such progressions may be 
based upon the relation q'^Xj = Aj+at giving Xpjy^k = A^g^^, k = 1,...,N, p E Z. The 
continuous spectrum may arise in the form of infinite number of permitted bands concentrating 
near the zero. 

For instance, the described situation may take place if q^u{qx) = u{x). A solution of this 
equation is given by u{x) = h{x)/x'^, where h{x) is a function periodic in the logarithmic scale 
h{qx) = h{x). For h = const, the wave functions iIj{x) are given in terms of the Bessel functions. 
For h < —1/4 there is the "fall onto the center" phenomenon and under certain conditions 
discrete spectrum consists of one bilateral geometric progression of negative eigenvalues p. 
The bad point of this and of the previous example is that the potentials u{x) are not bounded 
from below. Let us impose such a demand and assume that u{x) is bounded from below and 
non-singular. Evidently this condition does not permit the infinite negative spectrum, i.e. it is 
necessary to delete the lower tails of our geometric (or arithmetic) progressions. However, the 
formula q'^Xj = Aj+at still works and gives XpN+k = Xuq^'^ ,p = 0, 1, . . .. We have thus a partially 
self-similar spectrum: 

{g2A,} = {A,}/(Ai,...,A^), 

i.e. the scaling x ^ qx leads to the deletion of lowest eigenvalues from the spectrum. Below 
we outline a procedure for building a class of potentials with this property. 

Note that for the q^ ^ 1 case the geometric progressions of positive eigenvalues of the 
abstract operator L are bounded from below and no truncation is necessary. However, one has 
to go beyond the Schrodinger equation for building systems with such completely self-similar 
spectrum. 

3. Infinite-soliton systems. The inverse scattering method and the general theory of 
solitons (see, e.g., @, H, 0]) provide us a guide for construction of potentials with prefixed 
spectral properties. The compatibility condition of two linear equations 

Lip{x,t) = Xip{x,t), ipt{x,t) = Ailj{x,t) 

has the operator form Lf = [A, L]. When L = —d^ + u{x, t) and A = —d^ + 6ud^ + Su^ one 
arrives at the KdV equation ut + Uxxx — QuUx = 0. 

If the potential u{x) for x —>■ oo sufficiently fast then according to the inverse scattering 
method the potential u{x) can be uniquely reconstructed from the reflection coefficient R{k), 
discrete spectrum eigenvalues {A^} and the normalization constants {Of^} for corresponding 
eigenfunctions {ipjlx)}. For R{k) = and finite number of eigenvalues one arrives at the A^- 
soliton potentials. The corresponding solution of the KdV equation can be represented in the 
form 

u{x,t) = -2dl\nTNix,t), (2) 
where r^r is the determinant of a A^ x A^ matrix C, 

Tj, = det C, = S,, + ^^e(^»+^^)/^ (3) 

Ki ~r f^j 

Oi = hx - kft + 9f\ i, J = 1, 2, . . . , AT. 

Parameters ki describe amplitudes of sohtons related to eigenvalues of L as A,; 
are zero time phases of solitons and kf are their velocities. 



3 



Let us extend this solution to tlie pA^-solitonic one by imposing the constraint fc^+iv = 
qkj,q < 1, and take the p ^ oo hmit. This leads to the self-similar spectrum kpN+m = 
q^km, m = 1, . . . ,N, i.e. to the infinte-soliton wavetrains whose amplitudes form geometric 
progressions. The question of convergency of such a limit requires a special consideration. 

Applying scaling transformations to derived potentials one has to scale the time t too, which 
gives u{x, t) = q'^u{qx, qH). The phases of this potential are 6j{qx, qH) == kj+NX — kj^jsft + Of^. 
There are two "drawbacks" of the resulting expression. First, there are infinitely many arbitrary 
parameters 6f'\ i.e. there is a functional freedom in the definition of u{x). Second, there is 
a mixup between soliton amplitudes and phases — the solitons with the amplitudes kj^N 
have now the zero time phases of instead of e\%. Both these shortcomings are removed 

if we impose the constraint Of^N = Of'\ i.e. 6j{qx,qH) = 9j^N{x,t), which leads to truly 
self-similar infinite-soliton systems — the dilation x — > qx,t q^t just deletes solitons, 
corresponding to the lowest eigenvalues of L. Such potential is contains a finite number of 
parameters /ci, . . . , /cat, g, 9'i\ . . . , 9^^ and should be characterized by some "finite" equations. 
There are two equivalent approaches to the construction of these equations. The algebraic 
one, based upon the factorization method, and the functional-analytic one, relying upon the 
Darboux transformations. 

4. The factorization method. This method is well known in quantum mechanics 0, it 
was discussed already by Schrodinger. Within this approach one takes a chain of Schrodinger 
operators, Lj = —(P/dx'^ + Uj{x), and factorizes them as products of the first-order differential 
operators up to some constants Xf 

Lj = A+Aj + Xj, = T d/dx + fj{x). (4) 

Then the neighboring Lj are tied to each other through the abstract factorization chain 

Lj+i = A'^^-^Aj^-Y + Aj_|_i = Aj A'j + Xj, 

i.e. one passes from Lj to Ljj^i just by the permutation of operator factors. This gives the 
intertwining relations AjLj = Lj^iAj , LjA'j = A'j'Lj^i. As a result, if one has Ljip^^^ = 
Xip^^\ then ^(^+1) oc Aji^^l Indeed, Lj+iiAjip^^'^) = AjLjip^^'^ = XiAjip^^^). Under particular 
conditions the action of the operator Aj leads to removing or addition of an eigenvalue in the 
spectrum of the operator L^+i with respect to the spectrum of Lj. There may take place an 
isospectral situation as well. If zero modes of all Aj are normalizable, then Aj, Aj+i, . . . form 
the discrete spectrum of Lj. 
Let us define the operators 

Mj = Aj_^_j^_^ . . . Aj_^_^Aj , Mj~ = A'jA^^^...A'j_^]y_-^^. 

They relate eigenfunctions of the operators Lj and Lj^j^ due to the intertwining relations 

L,+nM- = M-L^, M+Lj+N = LjMf. 

It is easy to see that the products of operators Mf should commute either with Lj or with 
Lj+iy. Indeed, one has the equalities 

Af-1 Af-1 

M+M- = n (L, - A,+,), M-m; = n (W - (5) 

k=0 k=0 

which, together with the previous relations, look almost as an algebra of symmetries. 
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The main advantage of this method is that it is not tied to any particular type of spec- 
tral problem. One may realize the abstract factorization chain with the help of differential, 
difference, integral, etc operators. 

5. Darboux transformations. One starts from the semi-disrecte LA-pair 



= A-^)^^^ , A- = d/dx + fj (x) . 

The compatibility condition yields the intertwining relation AjLj = Lj^iAj , the resolution of 
which gives Uj = fJ — fjx + \j,Uj+i = Uj + 2fjx = ff + fjx + Xj- The substitution fj = — ^^r'Y^^-'^ 
converts the relation between Uj and fj into the equation + Uj(j)^^'^ = Xj(j)^^\ As a result 

one gets the original Darboux transformation: ■ip^^'^^^ = ip^'^ — {(p^^^ (f)^^^)-ip^^\ 

The compatibility condition (or the factorization chain) is equivalent to the following infinite 
chain of nonlinear differential-difference equations 

ifji^) + fj+iix))x + ffi^) - ff+ii^) = f^j = ^j+i - ^j- (6) 

One may search solutions of this equation in the form of power series in j. The finite term 
expansion occurs only if fj{x) = a{x)j+b{x)+c{x)/j, where a, b, c are some elementary functions 
of X 0. The resulting Schrodinger equation is solved in terms of the 2-^1 hypergeometric 
function. 

6. Polynomial supersymmetry. Upon the two neighboring operators Lj, Lj+i of the 
factorization chain a simple superalgebra describing the boson-fermion symmetry is realized. 
Let us introduce the 2x2 matrix Hamiltonian H = —d'^/dx'^ + ff{x) — fjx{x)(J'i and de- 
fine the supercharges = Af(ai ± 102) {ok are the Pauli matrices). Then one has 
{Q^,Q~} = H, [H,Q'^] = (Q^)^ = 0. This construction was generalized in to a sym- 
metry between particles with parastatistics (para-super symmetry). In the simplest case one 
takes the Hamiltoninan given by 3 x 3 diagonal matrix with the entries Lj^i, Lj, Lj+i. A 
symmetry between the boson and a parafermion of the second order is described now by some 
polynomial (cubic) relations between H and the correspodning generalized supercharges. 

As a next step in the "nonlinearization" of supersymmetry the following polynomial super- 
symmetry algebra has been derived in |^ 

{Q+, Q-} = PiH), [H, Q±] = (Q±)2 = 0, 

where Q± = Mtiai ± ia2)/2 and P{H) is a polynomial of H: 

N 

PiH) = l[iH-\k), H = 

k=l 

Since the supercharges have now more than one zero mode, the Witten index does not charac- 
terize the supersymmetry breaking. One may deform this algebra via the affine transformations 
and impose various natural constraints upon H. This gives an alternative way of derivation 
(which we are not describing here) of self-similar potentials considered in the next section. 

7. Self-similar potentials. The infinite soliton potentials with self-similar spectra consid- 
ered above appear as self-similar solutions of the chain (^. The simplest N = 1 case has been 
considered in |P, The general class of these potentials is defined by the following g-periodic 




reduction 111 



fj+N{x) = qfjiqx), iij+N = q l^j- (7) 
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Note that the sohtonic interpretation is vahd only if < < 1 and there are no singularities 
on the whole line —oo < x < oo, which is not true for arbitrary q,x,Xj G C. What are the 
properties of these potentials ? 

Algebraically, these are the systems whose symmetries are described by the quantum alge- 
bras. Indeed, the equations form a closed algebra in the cases when the operators Lj and 
Lj+7v are related to each other in some way, e.g., if they are similar to each other Lj ~ Lj^^. 
The g-periodic closure corresponds to the following constraint Lj+n = q'^TLjT~^, where T is 
the scaling operator TiIj{x) = \f\q\ip{qx). This gives Mj+Ar(x) = q^Uj{qx). 

For real g 7^ 0, T is the unitary operator T'^ = T~^. Substituting the operator constraint 
into (I) one arrives at the following polynomial algebra 

N N 

LB^ = q^^B^L, B+B- = X{{L-\), B' B+ = X{{q^ L - X^) , (8) 

k=l k=l 

where we have denoted L = Li, B^ = T^^Mi, B^ = M^T. For = 1 this is a g-analog of the 
harmonic oscillator algebra for = 2 one gets the SMq(l, 1) algebra, etc. 

Analytically, we deal with a class of functions appearing from solutions of the following 
system of nonlinear differential-g-difference equations: 

{fl{x) + f2{x))a: + /i (x) - = /ii, 

(/iv(a;) + qfiiqx))^ + f^{x) - q^ fl{qx) = /i^. 

A flavor of the structure of the general solution is obtained from consideration of various 
limiting cases, when fj{x) are expressed through known functions. The simplest situation is 
obtained when fj{x) = Cj = const, and arbitrary q. This leads to L = —dP/dx"^, i.e. to the 
free nonrelativistic particle, which is interpreted in this way as a g-algebraic system. When 
is not singular at a; = 0, the crystal base limit g results in the general A^-soliton 
potential. One can take the limit g — ^ 1 in such a way that the geometric progressions in the 
spectrum are converted into the arithmetic ones. Then for A^ = 1, 2 one gets the potentials 
u{x) oc x^, ax"^ + h/x^. For A^ = 3, 4 the functions fj{x) are expressed through the Painleve IV 
and V transcendents |Jl3|, [l^. For g = — 1 one gets a similar situation but now the functions 
are constrained to obey certain structure under the parity transformation. It is possible to 
consider the limit from Schrodinger equation to classical mechanics, in which case the self- 
similar potentials are determined by the much more simple functions |]TB|, |TB|. 



For g a primitive root of unity, g" = l,g 7^ ±1, one gets the finite gap potentials with 
additional (quasi) crystallographic symmetries ||T^. The appearance of Painleve functions and 



of the finite-gap potentials in a similar setting was discussed also in [|T8| , [T9[. It is natural to 
refer to the whole class of functions emerging for general initial conditions and < |g| < 1 
as the continuous g-Painleve functions. Because of the emergence of quantum algebras, the 
derived class of self-similar functions represents a new set of g-special functions defined upon 



the differential equations. They differ in structure from the basic hypergeometric series 
but there is a connection with them via the coherent states. 

Let us remark that the Darboux transformations allow one to remove or add only a finite 
number of levels from/to the spectrum of a given Schrodinger operator. If one finds an in- 
tertwining operator deleting or adding permitted bands from/to the given spectrum, then the 
systems with a self-similar infinite gap spectrum can be constructed from the requirement that 
the action of this intertwining operator is equivalent to a simple dilation of the coordinate x. 
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8. Coherent states. There are many definitions of coherent states in physics. The purely 
group-theoretical approach uses the orbits of groups lying behind the symmetries of a taken 
physical system. For simple spectrum generating algebras, as in the harmonic oscillator or 
singular oscillator cases, these states may be defined as eigenf unctions of the lowering operator. 
In this approach coherent states play the role of generating functions of irreducible represen- 
tations of the underlying algebra. Below we outline briefiy the results of application of the 
latter definition to the nonlinear algebra of symmetries (||). A more detailed discussion of the 



structure of these coherent states can be found in I^T . 

The action of operators upon abstract eigenstates of the operator L, L\X) = A|A), has 
the form: 

N N 

B-\\) = n VA-Afe|Ag-2), 5+IA) = Y[ ^Xq^-^Xq^). 

k=l k=l 

Let < g < 1 and A^ < 0. Then the A < eigenstates of L are formed from up to 
lowest weight unitary irreducible representations of the algebra 1^. Indeed, B~ is the lowering 
operator for A < states and a/A — A^ becomes complex for A < min{Afc}. Since B~\Xk) = 0, 
this problem does not arise for a special choice of A. Namely, if the point A < belongs 
to the spectrum of L, it must be of the form XpN+k = A^g^^ and the corresponding state is 
\^pN+k) oc {B^Y \ Xk) (we assume that |Afc) are normalizable) . 

Coherent states for these series are defined as eigenstates of the lowering operator B~: 

= a|a)W^ k = l,...,N. (9) 

There are N such states since they are defined for each lowest weight series separately. Repre- 
senting as a superposition of the states \XpN+k) one finds 



oo 




p=0 

where nVn-i is a basic hypergeometric series with the operator argument z = {—l)^aB^/ Xi . . . ) 
and the parameters {q'^,b'j} = {q'^Xk/Xj}. The general definition of the series of this type is 



|20 




bs J {q,bi,...,bs;q)n 

where r and s are arbitrary positive integers, and ai, . . . ,ar,bi, . . . ,bs are free parameters. We 
use also the following compact notations (a; g)„ = (1 — a) . . . (1 — and (ai, 02, ... , a^; q)n = 
(oi; q)n{ci2] <l)n ■ ■ ■ {o,n', q)n- These coherent states are normalizable if the complex variable a lies 
inside of the circle |ap < |Ai...AAr|. General ^V^s series are related to coherent states of the 
rational generalization of the algebra (|^) p2| . 

The A = eigenstates of L are simultaneously eigenstates of the B^ operators, i.e. they 
have a direct interpretation as coherent states. This is a special degenerate representation of 
the algebra (|). 

Unusual coherent states appear from the non-highest weight representations of (§) corre- 
sponding to the A > eigenstates of L. Let -Eq > be a discrete spectrum point. Then the 
action of B^ generates discrete spectrum of L in the form of one bilateral geometric progression 
Eoq"^"", n & Z (such situation cannot take place for the Schrodinger equation). In principle there 
may be an arbitrary number of such progressions because v^A — A^ is real for arbitrary A > 
and there are no truncation conditions as in the A < case. 
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Since for A > the lowering operator is B~^, the coherent states should be defined as 
eigenstates of this operator (instead of B~): 

B+\a)+ = a\a) + . (10) 

These coherent states can be expanded over the series of eigenstates of L with positive eigen- 
values |£'oQ'^"): 

°° / 0, . . . , \ 

oc Cna"\Eoq^") Qc oi^N [ . , ; q'^ , z ] \Eo) , 



^^i/Eo, . . . , Aat/ Eq 

where o'^Pn is a bilateral g-hypergeometric series with the operator argument z = aB^ / (— -Ec 



The general representative of these series is defined as follows [|20 



ybi,...,bs J ^t^^ {bi,...,b,;q)n ^ ^ 



Coherent states are normalizable if the parameter a lies outside of the previous indicated 
circle |ap > |Ai . . . AtvI- 

When the A > region is occupied by continuous spectrum, coherent states are defined as 
eigenfunctions of B~^ again, but in this case it is necessary to use integrals in the expansion of 
|«)+ over the states |A). Since the continuous spectrum may be considered as a continuous direct 
sum of irreducible representations with bilateral geometric progressions formed by spectral 
points, there are now infinitely many coherent states. Let = 1 and the states |A) are not 
degenerate. Then, normalizing Ai = l/(g^ — 1), we have 

\a)\ = C [a) ' 



-Ag2(l - q^)-q 



2\ 

oo 



where 

27iis — ln(ag^-\/l — q'^) 
In 

These states are normalizable for |ap > 1/(1 — g^), and have the unit norm for 



7. = r-^ , seZ. 



(-Ag2(l - g2).g2)^ sinvrr iq'^;q'^] 



oo 



where r = In \aqy/l — g^l/lng. The last integral is calculated exactly being a particular subcase 



of a Ramanujan g-beta integral |20 



Schrodinger operators with self-similar potentials have the continuous spectrum for A > 0. It 
is doubly degenerate, which leads to duplication of the number of coherent states. An instructive 
example is provided by the free nonrelativistic particle model for which L = —(P/dx^. The 
corresponding g-harmonic oscillator algebra generators have the form 



B- = T-\d/dx + 1/^1 - g2), B+ = {-d/dx + l/\/l - g2)T, 

B-B+ -q^B+B- = 1, L = B+B- - l/{l-q^). 

The eigenvalue problem B''tp'^{x) = a4'^{x) leads to the pantograph equation, which was 
analyzed in detail in |23]. The corresponding results show that in accordance with the purely 
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algebraic consideration there are no normalizable coherent states of this type. However, there 
are infinitely many normalizable eigenstates of the operator determined as appropriate 
solutions of the advanced pantograph equation 



An important property of these states is that they are determined by C°° but not analytical 
functions at the x = point. Consideration of the > 1 symmetry algebras in this realization 



leads to the generalized pantograph equations p4| 



As a general conclusion to this section we would like to stress that all three types of objects 



associated with g-special functions in their standard meaning the ordinary and bilateral 
basic hypergeometric series, and Ramanujan type integrals, show up in the context of self- 
similar spectral problems. 

9. Second order difference equation. Let us describe briefiy a realization of the algebra 
of symmetries (H) upon the discrete Schrodinger equation or three term recurrence relation for 
orthogonal polynomials [^, One considers an infinite chain of Jacobi matrices Lj and 

the corresponding eigenvalue problems: 

Lj^i = + + Ki^i = ^^L e Z. (11) 

In the case of orthogonal polynomials, one considers ( ]TT[ ) only for n > and imposes the 
boundary conditions ipoiX) = 0, ipKX) = \ — b{. Simlar to the continuous case the forward 
discrete time step is defined by the Christoffel's transformation to kernel polynomials 

A — Xj+i A — Xj+i 
The backward transformation was analyzed by Geronimus and it has the form: 



^Pt' = iji + AiijL, ^ R,^i. (13) 

Here A{ and Cl are discrete analogs of the superpotentials fj{x) and Sj, Rj are the first order 
difference operators. The compatibility conditions of the j — j±l moves yield the factorizations 
Lj = S jRj -\- Xj ^ Lj—i = RjS j -\- Xj ^ OTui^ = Ai^Cl, bl = A-l^_^^+Cl+Xj. The abstract factorization 
chain in this case is equivavlent to the following set of nonlinear finite-difference equations: 

Mi^n-l = A^n ^^i, "^5 A^n~^ Cn~^ 'Xj = ^n+1 ~^ ^ + Aj_i, (14) 

known to define a discrete-time Toda lattice. Note that this is not an isospectral fiow since the 
constants Aj determine the character of the change of the spectrum of the taken Jacobi matrix 
after the Christoffel or Geronimus transformations. 

Discrete Schrodinger equation analogs of the self-similar reductions (|^) were described in 

= = qCU, A,+^ = gA„ (15) 

where k is an integer. When n is considered as continuous, k may be taken as a continu- 
ous variable as well. This closure is associated with classical, semi-classical and, so-called. 



Laguerre-Hahn class of orthogonal polynomials on linear and g-linear grids The corre- 



sponding recurrence coefficients are related to ordinary and g-analogs of some discrete Painleve 



transcendents |2^, The algebra of symmetries of these systems is derived along the same 
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lines as in the continuous case, however, in the present case the compact versions of algebras, 
like suq{2), are allowed as well. The coherent states are defined in the same way as in the 
continuous case. For some explicit examples see, e.g., P2 . 



The following discrete-time Volterra lattice has been derived in 



Di 



-1 



(16) 



It can be mapped upon the discrete-time Toda lattice (jT^) via the following quadratic relation: 

Xj = const. — Pj, 

which generalizes well-known relation between the ordinary Toda and Volterra lattices. There 
is also the second similar mapping 



Ai 



Cri — {.D2n — l3j){D'2n+l ~ Pj), 



(18) 



with the same connection between Xj and Pj as given before. As shown in this discrete- 
time Volterra lattice is related to the (^-algorithm proposed by Bauer in numerical analysis 
33| . There exists an ansatz of semi-separation of discrete variables in (0) which leads p6| 
to recurrence coefficients of the Askey- Wilson polynomials — the most general set of classical 
orthogonal polynomials |Q and, simultaneously, of the polynomials considered by Askey and 
Ismail . 

An interesting discrete symmetry for the chain 



has been described in 

with a freedom in the intermediate steps of two-step discrete-time shifts j 



It is associated 
j + 2 and related 



to the statement on permutability of a sequence of Backlund-Darboux transformations for fixed 



set of corresponding parameters. In |^ an analog of this symmetry for the discrete time Toda 
and Volterra lattices has been derived. 

Let us describe briefiy this refactorization symmetry for the equation (|16|). Let jSj ^ for 
any j, then: 



f3,Di + 



n+l 



n-l 



^^+i)+/9,-i(/5,- 



D 




(19) 



n 



-1 



n-2 



Di 



/3,-i(/?, 



The change of spectral parameters (3j looks as follows: 



(20) 



(21) 



If one substitutes into (|T6|) instead of Df^, Dl~^ and the tilded variables and keeps all 

other D'^, P^, k ^ j, j — 1, fixed, then the resulting equation will be satisfied automatically, i.e. 
we have a discrete symmetry. The transformation laws for superpotentials y4{, C^, describing a 
similar symmetry for the discrete-time Toda lattice, follow from the maps (|l^, (piSl) . If Pj = 
or Pj-i = (j is fixed), then there appears some additional freedom, we refer for details to 

We have sketched only the simplest types of discrete-time integrable systems. For an analysis 
of ordinary Toda lattice and its self-similar reductions see, e.g. 
examples and their applications see [I 



36| , [37[] . For more complicated 



39|, HOl El|. Recently self-similar reductions of spectral 



transformation chains have been considered for biorthogonal rational functions in |^ . 
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10. One-dimensional Ising chains. Tau-function of the iV-sohton solution of the KdV 
equation (^) can be represented in the following Hirota form, which was widely discussed in 
the literature (see, e.g., |^3|, Q and references therein): 

( ^ \ 

tn = exp \ Aijfiifij + OiHi . (22) 

tii=0,l \l<i<j<N i=l J 

Here the phase shifts Aij are expressed via the spectral variables ki in a simple way 

^A., ^ (^iJlM!. (23) 



{ki -\- kj 



In it was noticed that this expression for r^v defines the grand partition function of the 
lattice gas model if 9i = 9^^^ = const. A comparison with shows that /Xj play the role of the 
filling factors of the lattice sites by molecules and 9^^^ is a chemical potential. The constants 
Aij are proportional to the interaction energies of molecules. 

Simultaneously, tn describes the partition function of a particular one-dimensional Ising 
chain: 

Zn= Y: e~^^ (3=^, (24) 

tTi=±l 

N 

E= Y Jij(^i(^j -Y^i^i^ 

l<i<j<N 1=1 

where N is the number of spins cTj = ±1, Jij are the exchange constants. Hi is an external 
magnetic field, T is the temperature and k is the Boltzmann constant. Indeed, after the 



substitution into ([22|) of the spin variables /ij = (o-j + l)/2, some simple calculations yield 

1 1 ^ 

with the following relations between soliton parameters and Ising chain characteristics 

1 11^ 
/3J,, = --A,„ (3H, = -9, + - Y A,,. (26) 

A similar relation with Ising chains is valid for the whole KP hierarchy and some other 
differential and difference nonlinear integrable equations. The corresponding tau-functions 
have the form (^2]) for different choice of the phase shifts Aij and the phases 6'^, a list of such 
equations can be found, e.g., in |4^, ^. 



A crucial observation of |^4[ is that the condition of self-similarity of the spectrum is related 
to translational invariance of the spin exchange constants Jij of the infinite Ising chain induced 
by the KdV equation. In the simplest case one demands translational invariance of the system 
with respect to the shift by one site j — > j + 1 which means that Jj+ij+i = Jij. As a result, 
the exchange Jij or phase shifts Aij depend only on the distance between the sites |« — j|. This 
natural physical requirement forces ki to form one geometric progression 

ki = kiq'-\ q = e-^'', = 2 In | tanha(z - j)|, (27) 

where ki and g < 1 are free parameters. One may demand also the invariance with respect 
to shifts by a multiple of the lattice site, when Jj+mj+m = Jij, and this leads to the general 
self-similar spectra kj+M = qkj. 
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In fact the translational invariance in not exact for finite chain 1 < j < N. The hmit 
N ^ oo corresponds to the thermodynamic hmit. It gives an infinite soliton potential with self- 
similar discrete spectra. The coordinate x and time t (and higher "times" of the corresponding 
hierarchies) are interpreted as parameters of the magnetic field H^. The x,t dependent part of 
Hi decays exponentially fast for i oo, since q < 1. Therefore in the thermodynamic limit 
only the values of constants Of'^ are relevant for the partition function. The formalism allows 
us to treat the M-periodic magnetic fields, -ffj+M = Hi, which for M = 1 is just a homogeinity 
condition. 

For the KdV equation case one has < |tanha(i — j)| < 1 and = —Aij/ip > 0, i.e. 
an antiferromagnetic Ising chain. A similar situation takes place for general M-periodic case. 
Such interaction has a long distance character but the intensity falls off exponentially fast. The 
absence of phase transitions in such systems for nonzero temperatures is well known. 

For a — i> or g — i> 1 the phse shifts Aij oc Jij/kT are diverging. If one renormalizes the 
exchange constants J^^^"^ = Jij{q~^ — q) and the temperature kT^en = kT{q~^ — q), then the 
interaction energy of a single spin with all other ones will be finite for g — 1. As a result, in 
this limit one actually gets a nonlocal interaction model with a low effective temperature. Note 
that for imitation of the change of the temperature it is necessary to change simultaneously the 
magnetic field H = if""^"/ — g). 

The limit g ^ gives J^^"' oc or the high temperature nearest neighbor interaction 



Ising chain (if H is renormalized) . If H is kept finite then this limit corresponds to the non- 
interacting spins. So, the formalism provides only a partial description of the partition function 
in a two-dimensional subspace of (T, H, q) . For fixed q the temperature T has a prescribed value 
and one may normalize the "KdV temperature" to /3 = 1. 

Using the Wronskian formula for the representation of A^-soliton potential , we were able 



to calculate the partition function for the translationally invariant Ising chains in a homogeneous 



magnetic field. Omitting technical details, which can be found in |^^, for M = 1 we get 



Zn exp(— A^/?//) for A^ oo, where the free energy per site // has the form 

- (3fj{q, H) = In + ^ r duH\piu)\' - gtanh^ (3H), {2i 

(g2; q^vl^ 47r Jo 



/oo 



'^^ ^' (gV-g^)L(g%-^'';g^)L Asm\v/2) ^9l{v/2,q^y 

Here ^i,4(?/, q) are the standard Jacobi theta-functions of the argument y and base q (our base 
is (f). The original expression for the density function p(z/) is 

where a small regularization parameter e — > is introduced for absolute convergency of the 
infinite sum. The last equality in ( pQ]) was obtained with the help of the iipi Ramanujan's 



bilateral basic hypergeo metric series sum [20 



g («; <i)n _ (9> <l/<^Z- q)c 



{b;q)n {b,q/a,z,b/az;q)c<, ' 

This fact demonstrates again the general statement given earlier that the g-special functions 
are always behind the self-similar spectral problems. 



12 



The derivative of fi{H) with respect to H yields tlie total magnetization of the lattice: 



tanhi 



1 ^ 

m{H) = -dnfi = lim 

N^oo iV — ; 

4 = 1 

V ^|(z/,g2)cosh2/?if-^2^z/,g2)sinh2/?iyy ■ ^ ' 

The function m{H) grows monotonically with iJ, its quahtative properties are predicted by the 
general theorems on the behavior of ID systems with the fast decaying interactions. 

The main drawback of the described scheme is that the exact calculation is valid only for 
a fixed temperature in the partition function (g is fixed). One may try to replace ( P7D by 
Aij = 2n\n \ {ki — kj)/{ki + kj)\, where n is a sequence of integers, to change the temperature 
and look for integrable models with such phase shifts. This is not a simple task since there are 
many constraints involved into the resolution of this problem. We have found only one more 
example for n = 2 corresponding to a subcase of the iV-soliton solution of the KP-equation 
of B-type (BKP). The BKP equation has the form (see, e.g., 0, ||]) 



9lJ--5ir-ij-2+lJ-5 -30Tr-ij- + 30— _ 5— = 0. (31) 



dxi \ dx5 dx^dxi dxf dx^dxi dxidxf \dxi J j dx 

r-function of this integrable equation generates partition functions of quite general Ising chains. 
The exchange constants have now the form 

^ ^ (a, + a,)(6, + 6,)(a, + 6,)(6, + a,) ^""^^ 

with the energy and partition function given by (p^ . For Oj = hi = ki/2 one gets the KdV- 
inspired model for j3 = n = 2. 

Translational invariance of the spin lattice, Jij = J{i — j), results in the following spectral 
self-similarity 

ai = q'-\ bi = bq'-\ g = e'^", (33) 
where we set ai = 1 and assume that a > 0. This gives the exchange 



4 coth'a{i-j)-{b-iy/{b+iy 

Since this expression is invariant with respect to the inversion b b^^, the parameter b is 
restricted to the unit disk |6| < 1. It is not difficult to see that for —l<b<—q one has now 
the ferromagnet, i.e. Jij < 0. For two other physical regions q < b < 1 and 6 = e**^ 7^ —1, one 
has the antiferromagnet, i.e. Jij > 0. 

In the thermodynamic limit — > 00 the partition function can be represented as the 
determinant of a Toeplitz matrix, which is diagonalized by the discrete Fourier transformation. 
Using again the iipi sum, we have found the free energy per site for the homogeneous magnetic 
field Hi = H: 

-Pfi[H) = -\n- +1- duln\2p{u)\, 

4 {-q,-q,-bq,-q b;q}^ An Jo 
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where 



(-e^", -ge— ; g)oo V {h-\qh-q)^ {b,qb-\q)^ J 

The magnetization if obtained after taking the derivative with respect to H 

1 - - / ^-777 , (35) 

vr io 1 + d{u) cosh 2f3H / ' ^ ^ 

where 



for g < & < 1, 

for — 1 < 6 < — g, and 

d{u) = - 



{-q; g)^2Im Oi {v - {i/2) In 6, g V2) 

" (-g;g)^2Re^2(i^-(V2)ln|6|,gi/2) 
(ge^^ ge-^'^; g)oo2 sin(0/2) g2(^^, q^'^) 

-q; q)loM'^ + 0/2, - - 0/2, gV^)] 



for b = e**^. Here ^2(i^, Q^^"^) is another Jacobi 6'-function. Taking the hmit 6 — 1 one may find 
the magnetization for the second value of discrete temperature of the KdV equation inspired 
spin chain, n = 2. The details of calculations and graphical representation of m{H) for some 



values of the parameters can be found in ^J] . In a somewhat similar way one may consider the 
M > 1 periodic systems as well. 

11. Random matrices. Random matrices were employed by Wigner and Dyson for 
studying spectra of complex systems with many degrees of freedom. When the density of 
levels is high enough, excitation energies can be described statistically. The detailed structure 
of hamiltonians is not known and statistical characteristics of the system are described by 
averaging over ensembles of random matrices. A fundamental constraint is that the probability 
distributions should be invariant under basic symmetries such as parity, rotation and time- 
reversal transformations E9| . 



Gaussian ensembles use a real symmetric, Hermitian or self-dual Hermitian random matrix 
H with the statistically independent elements ifj^. The probability density P{H) for Hi^ to lie 
in a unit volume is proportional to exp {—atiH^ + bti H) , with the measure invariant under 
orthogonal, unitary or symplectic transformations. Abandoning the condition of statistical 
independence one may pass to ensembles with more complicated P{H). 

Dyson has introduced circular ensembles of unitary random n x n matrices 5* with eigen- 
values ej = e*'^^, j = 1, . . . ,n, such that the behavior of the phases 0j is equivalent to the 
distribution of eigenvalues of a system. The condition of invariance of the measure under all 
unitary authomorphisms S USW, where U, W are arbitrary unitary matrices, determines 
the probability distribution. Any S can be represented in the form S = U~^EU, E a diagonal 
matrix with the eigenvalues ej and U some unitary matrix. Then, using the invariance of prob- 
ability measure, one can gauge away the [/-dependent part of S and integrate out the "angular 
variables" in the probability distribution. As a result, the eigenvalue distribution becomes 



Pdcpi . . . dcpn Y[ 

i<j 
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One can weaken the condition of invariance of the measure under arbitrary left or right transla- 
tions and leave only pure unitary transformations. Then the invariant measure is not uniquely 
defined and the eigenvalue distribution may be written as 



/(0i,...,0„)n 

i<j 



Pdcpi . . . d(j). 

where the function / is symmetric under permutation of its arguments. One of such ensem- 
bles has been considered by Gaudin long time ago |50|. In this model one has the following 
probability law for eigenvalues 



Pd(t)i ... I 



Kj 



€2 6 J 



(36) 



interpolating between the Dyson unitary ensemble {z = 0) and the uniform distribution {z = 1). 

The same model can be interpreted also as a Coulomb gas on a circle with the partition 
function 



y f f 

Zn OC ... I 

JO Jo 

where (3 = 1/ kT is fixed and the potential energy is 

sinh^ 7 



exp(-/5^r(0, 

i<j 



(3V{(j)i - 0,) = In 1 + 



sin^((0,-0,)/2). 



-27 



(37) 



The grand partition function of this model corresponds to the r-function of the KP hierarchy 
in the specific infinite soliton limit |4^|. The finite soliton solutions provide thus a discretization 
of the system, namely, a lattice gas on the circle. For example, one may take unitary n x n 
matrices S with the eigenvalues equal to iV-th roots of unity, i.e. = 2Trmj/N, rrij = 1, . . . , N. 
The measure is taken to be continuous in the "angular" variables and discrete in eigenvalues, 
i.e. one takes the sums over (pj instead of the integrals. The completely continuous model is 
recovered for N ^ 00, (pj fixed: 



n N 



N 



(I) 



N^co JO 



mi=l m„=l 

For n X n matrices the lattice partition function is 



2n 



2n 



(38) 



Zn{N,z) 



27r 



N N 

EE n 

mi=l m„=l l<i<j<n 



The grand canonical ensemble partition function is 



2i'Kmj 



-27 



Ziz, 



E 



N 



nl 



^ exp J2 ^mkfJ'ml^k + {9 + V) J2 



(39) 



m=l 



n=0 Atm=0,l \l<m<fc<Af 

where rj = ln(27r/A^) enters as an addition to the chemical potential 6 and 

sin^(7r(m - A;)/A^) {am - ak) {hm - bk 



A. 



mk 



In 



In- 



sin^(7r(m - k)/N) + sinh^7 (a^ + 6fe)(&m + Ofc) ' 
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is the KP phase shift with the restricted choice of parameters 

a^ = e2'™/^, bm = -zam, m = 1, 2, . . . , iV. (40) 

So, the grand partition function of this discrete circular unitary matrix model coincides with 
the particular KP A^-soliton r-function at zero hierarchy "times". In the thermodynamical 
limit N —>■ oo the relation (^) takes place and we get the matrix model (pG]). In December 
1998 we have known that this root-of-unity discretization was considered by Gaudin himself 



51| , where the connection with Ising chains was noticed as well but the relation with integrable 
equations was not established. 

The BKP equation suggests a further generalization of this matrix model with the proba- 
bility law 



Pdcpi . . . dcpn oc 



+ 9 



ti + Ztj 



(41) 



The grand partition function is also defined by (^) where A^k is given by (0). In order to 
escape singularities appearing from zeros + e^- = for (pj = 0j + tt, the parameters am,bm 
have to be restricted. The choice (^) works for odd N only, which makes the N —>■ oo limit 
problematic. Another option is to replace m/N in (^) by m/2N, which corresponds to the 
half-circle system. Comparing the connection of BKP equation with Ising chains described in 
the previous section we see, that now we have spins on the circle, g is a root of unity and the 
exchange Jmk = Jijn — k) is 

^ ^ , . 11 + sinh^ 7/ cos^inm/N) 

pj[m) = —-In -2 — ; 27 ; — ^• 

4 1 + sinh 7/ sin [irm/N) 

An important difference between the < 1 and |g| = 1 Ising chains is that in the former case 
the thermodynamic limit N ^ 00 does not lead necessarily to the continuous space models as in 
the circular case. Also the relation of non-compact Ising chains to hermitian random matrices 
is not so straightforward. In particular, in order to get nontrivial weight functions one has to 
take inhomogeneous magnetic field corresponding to nonzero hierarchy "times" . Calculation of 
the partition function is more involved in this situation. 

12. Conclusions. Self-similarity in spectral problems serves as a key to special functions, 
in particular, to g-special functions. We hope that the brief review given in this paper gives 
a sufficient number of arguments for convincing in this vague statement. As to the notion of 
exact solvability of a given equation, it requires a more rigorous group-theoretical definition in 
terms of the operator (differential, difference, etc) Galois theory. In this approach one fixes in 
advance the basic fields of functions entering spectral problems and asks about the transcen- 
dence of solutions of the resulting equations over these fields. In the scheme we have presented 
the situation is slightly different — the basic fields of functions are defined as solutions of com- 
plicated nonlinear equations appearing after self-similar reductions of some infinite-dimensional 
systems of equations. Usually this gives simultaneously a part of solutions of the underlying 
linear spectral equations in quadratures over the basic field of functions, but determination of 
the trancendence of the general solution remains open. Therefore these two approaches should 
be used in conjunction. One of the problems appearing on this route is to give a complete char- 
acterization of the changes in the structure of corresponding (differential, . . . ) Galois groups 
resulting from the Darboux-type spectral transformations ||21|| . 

Self-similar spectra considered here are based upon linear map of the spectral parameter 
X ^ qX + h. This leads naturally to spectra containing a finite number of geometric progressions 
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and "Painleve-type" transcendental functions on the background. A natural generalization in- 
volves polynomial maps A -P(A) leading in general to chaotic structure of the spectrum. The 
situation now is much more complicated and only the simplest possibilities were considered in 
the literature. Note that, similar to the linear case, superpositions of several spectral sequences 
generated by the map A —>■ -P(A) are permitted too, but strong difficulties arise in the charac- 



terization of functions involed in these systems already in the simplest cases (see, e.g., [0, |53|). 
The scalings of discrete grids, or decimations entering at this stage hint upon the relevance of 
the self-similarity hidden in wavelets and quasicrystals. Perhaps all these subjects are naturally 
unified within an extended appoach to the concept of self-similarity. 
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